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TECHNICAL PAPER-1389 


SUMMARY OF ATMOSPHERIC WIND DESIGN CRITERIA FOR WIND 
ENERGY CONVERSION SYSTEM DEVELOPMENT 

Chapter 1.0 Introduction 

A highly condensed version of Chapters 2 through 5 on wind characteris- 
tics from the "Engineering Handbook on the Atmospheric Environmental Guide- 
lines for Use in Wind Turbine Generator Development, " NASA TP-1359, is 
presented in this report. Basic design values of the most significant wind 
criteria are presented in graphical format. The design values are given without 
discussion of the physical processes involved or of the analytical methods used 
to develop the design curves. For these details the reader should consult the 
previously mentioned engineering handbook. 



Chapter 2.0 Wind Speed 


2. 1 Extreme Wind Speed 

2. 1. 1 Extreme Wind Speed at a Height of 10 m 

The extreme wind speed at a height of 10 m, W (W, nn £l _) , is 

* h = 10m v h=30 ft' 

given by Figure 2-1 in terms of risk of exceedance. The engineer must select 
the degree of risk he is willing to accept that the extreme wind speed designed 
for will be exceeded at least once during the expected life of a Wind Turbine 
Generator (WTG) . The degree of risk conventionally used ranges from aero- 
space values of 10 percent for an expected life of 25 years to building code 
values of 63 percent for an expected life of 50 years. 

2. 1. 2 Adjustment of Extreme Wind Speed for Height 

The value of wind speed selected from Figure 2-1 is adjusted for height, 
h, by the adjustment factor, k^ (kp, given in Figure 2-2, i.e. , 


( 2 . 1 ) 

where 

W^ = extreme wind speed at height h 
k^ = adjustment factor for height h. 


W, = k W. 


h 


1 h=10 m 


< k i*h=30£t> 


2. 1. 3 Adjustment for Time of Structural Response to Gust 

The extreme wind speed, adjusted for height as described in Section 
2. 1.2, is further adjusted for the response time of the structure by multiplying 
by the factor k (k’ ) given in Figure 2-3. The adjustment factor is a function 

of the extreme wind at height h as determined from Sections 2.1.1 and 2. 1. 2. 


Three categories based on structure size are specified in this regard: 
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Category a — Structures or structural components of 20 m (65 ft) or 
less in extent in any dimension. 

Category b — Structures or structural components larger than 20 m 
( 65 ft) but for which neither the greatest lateral nor 
vertical dimension exceeds 50 m (165 ft). 

Category c — All structures larger than those in Category b; for 
Category c, k^= 1.0 (k^ = 1.0). 


2. 2 Mean Wind Speed 

2. 2. 1 Annual Mean Wind Speed 


The approximate areal distribution over the contiguous United States of 

annual mean wind speed equal to or greater than W. (W. _,) is given 

h=10 m ' h=30 ft' 

in Figure 2-4. The designer should select the design value of annual mean wind 
speed from this curve based on the percentage of the area of the country to 
which he anticipates sale of the WTG. 

2. 2. 2 Adjustment of Annual Mean Wind Speed with Height 


The adjustment of annual mean wind speed with height is achieved by 
multiplying the wind speed determined in Section 2. 2. 1 with the adjustment 
factor k (k') given in Figure 2-5, i.e. , 

O «J 



WlO m 


( k 3 W h=30 ft) 


( 2 . 2 ) 


where 

A 

W. = annual mean wind speed at height h 
h 

k = adjustment factor for height h. 


2. 2. 3 Wind Speed Duration Curve 

The number of hours, T, for which the wind speed, W^, is expected to be 
equal to or exceed a prescribed value, W , is estimated by 
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(2.3) 


A 2 

= 8766 exp [-7r(W p /2W h ) ] 


Any units can be employed for the wind speed provided they are consistent for 

W and W , respectively. A plot of the wind speed duration curve is given in 
p n 

Figure 2-6. 
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ADJUSTMENT FACTOR (DIMENSIONLESS) 


HEIGHT, h (m) 


2 5 10 50 100 500 800 



Figure 2-2. Factor for adjusting extreme wind speed with height. 
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ADJUSTMENT FACTOR (DIMENSIONLESS) 


W h (m s' 1 ) 



10 


20 


50 100 200 

W h (mph) 

Figure 2-3. Adjustment factor for response time of structure. 
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ADJUSTMENT FACTOR (DIMENSIONLESS) 
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Figure 2-5. Factor for adjusting annual mean wind speed with height. 
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Figure 2-6. Number of hours per year, T, wind speed is expected to be greater than 

Wp (wind speed duration curve) . 


Chapter 3.0 Mean Wind Speed Vertical Gradient 


The variation of horizontal, mean wind speed with height, AW, /Ah, is 

h 

presented in Figure 3-1. The wind shear is expressed in dimensionless form 
(h/W^) AW^/Ah and is given as a function of height made dimensionless with 

surface roughness length, z . Typical values of surface roughness length are 
given in Table 3-1. 


TABLE 3-1. TYPICAL VALUES OF SURFACE ROUGHNESS LENGTH 



Surface Roughness Length, z 

o 


(m) 

(ft) 

Sea or large bodies of water 

10 -5 - 10“ 4 

3 x 10 -5 - 3 x 10 -4 

Open country with no obstructions 

10 -4 - 10“ 2 

3 x 10 -4 - 3 x 10“ 2 

Open country with scattered windbreaks 

10 -2 - 10 -1 

3 x 10 -2 - 3 x HT 1 

Country with many windbreaks, small 
towns, outskirts of large cities 

ur 1 - i 

3 x 10 -1 - 3 

Surfaces with large and frequent 
obstructions (e.g. , city centers) 

1-4 

3 ~ 13 
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Figure 3-1. Dimensionless vertical gradient in horizontal mean wind speed. 


Chapter 4.0 Turbulence 


4. 1 Spectral Model 

4.1.1 Spectra 

The turbulence kinetic energy spectral densities for the atmospheric 
boundary layer (to elevations of 150 m) recommended for WTG design are-. 

Longitudinal 


0 



12m 3 W h=l0 m h[in ( 10/z o + ^ in ( h/z D + 1 )J~ 

1 + 192 [hftin(l0/z + l)/W_ in(h/z + 1) ] 5/3 

o ' h=l0 m v o ' 


Lateral 


0 



4 ‘ 0W h=10 m h t jen ( 10 ' /z o + 1) fn(h/Z Q + 1)] 

1 + 70 [hn*n(l0/z +1)/W, _ in(h/z +l)] 5/3 
o ' h=l0 m v o ' 


Vertical 


0-5W m h[in( 10/z + 1) in(h/z +1)] 

* (ft) = ^m ^ o ' ^ o } ± 

z 1 + 8[hn!n(lO/z Q + 1) / W h=1Q m in(h/z Q + 1) ] 7 


(4.1) 


where 

0 W ( n ) = spectral density distribution of turbulence kinetic energy 
a 

a = designates either x, y, or z component of the fluctuation w 
ft = frequency in cycles per second, Hz 


i&> 
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W, „ . = horizontal mean wind speed at h = 10 m 

b=10 m 

h = height above natural grade 

z = surface roughness length (typical values of surface roughness 
length, z q , are given in Table 3-1) . 

The terms h, z , and W. are expressed in meters and meters per second, 

o h=10 m 

respectively. The units for cp are then meters square per second. The sub- 

a 

scripts w , w , and w refer to the longitudinal, lateral, and vertical wind speed, 

x y z 

W . The coordinate system (x, y, z) is chosen with the x-axis oriented along 

the mean wind direction which is assumed to lie in a plane parallel to the Earth’s 
surface. 

A plot of the turbulence kinetic energy spectra in dimensionless coordi- 
nates is given in Figure 4-1. In dimensionless coordinates, the spectra for all 
three wind speed fluctuations lie on the same curve. 

4. 1.2 Turbulence Intensity (rms Value) 

The recommended value of turbulence intensity, a , is given by 

a 


<T = K W T, 
w 4 h=10 m 

a a 


where is given in Figure 4-2 as a function of surface roughness, z q , and 
a 

height, h. 

4. 2 Discrete Gust Model 
4. 2. 1 Extreme Discrete Gust 

4. 2. 1. 1 Shape of Longitudinal Extreme Discrete Gust 

The extreme, longitudinal discrete gust shape recommended for design 
is given by 
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1 /3 

W G U) = 1. 8 W.{l - exp[-(sin(7r£ > /2a)) ]} ; 0 < C, < a (4.2) 

1 /Q 

w _(£) = 1.8 W.{1 - exp[-(sin(7r(l - ^)/2(l - a))) ]} ; a < t, < 1. 0 

Ijr 1 


where 


£ = t/r ; a = 0. 12 + 0. 05 in h 

and 

W (£,) = instantaneous horizontal wind speed fluctuation about the 
mean, i.e., W Q (£) = W(C) - 

W. = the effective average horizontal, discrete gust amplitude, i.e. , 

1 w. = W(t) -W, 
i ' h 

W(t) = the average over the period r of the wind speed centered around 
the peak wind speed 

t = time 

r = period of gust 
h = height 

a = effective gust rise time. 


4. 2. 1.2 Gust Period 

The recommended period of the gust in seconds, r, is given in Figures 
4-3 and 4-4. The period t is expressed as a function of the spatial distance, 
Aa, over which the gust is 50-percent coherent (i. e. , the spatial distance over 
which the gust is estimated to be effectively uniform or the distance it would 
engulf) . Figure 4-3 gives the period of a gust effectively coherent over a hori- 
zontal distance, Aa = Ax, and Figure 4-4 gives the period of a gust effectively 
coherent over either the lateral or the vertical distance, Aa = Ay = Az, 
whichever is of interest. 
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4„2. 1.3 Extreme, Longitudinal Discrete Gust Amplitude, W 

Values of W are computed from \V. = where is given in Figures 

4-5 through 4-8 as a function of height, h, period, r, and mean wind speed at the 

10 -m level, W, . . The value W. represents the three standard deviation 

h=l0 m l 

value (99 percentile). That is, there is only an approximate 1-percent chance 
that W. will be greater than the value given in Figures 4-5 through 4-8. The 

value of is computed from W^_ ^ ^ given on the respective figure by the 

relationship 



k W, 

6 h=10 m 


(4.3) 


where k , which is a function of z and h, is determined from Figure 4-9. Thus 


W. 

l 


k k w 

5 6 


h= 10 m * 


4.2.2 Cyclic, Discrete Gust Model 
4. 2. 2.1 Cyclic, Discrete Gust Shape 

The conventional shape used for a cyclic, discrete gust is: 

(4,WJ = W. ( 1 — cos 2 7rf, ) ; 0<£<1 (4.4) 

Go h Ao? 


where 



and 


r = period of gust 

W = gust amplitude relative to the mean wind speed in the direction a 
A a 

W = instantaneous wind speed relative to the_mean wind speed in the 
" direction, ty ; i.e. , W^U) = W x (&) - W h> W^U) = W^), 

^d w Gy ( ^) = W y (0- 

In applying Equation 4.4 the value of t is determined from either Figure 4-3 or 

4-4. 

4. 2. 2. 2 Cyclic, Discrete Gust Amplitude 

The recommended discrete gust amplitude is 


W A = vk kk(W , n ) . 

An 4 a 6 7' h= 10 m ' 


This value of represents the gust amplitude average over a year based on 

a Rayleigh distribution of wind speed, and v is the ratio W. /a The value of 

A a ext 

v is selected by the designer and effectively corresponds to the number of stand- 
ard deviations represented by the imposed gust amplitude . The distribution 

of W, is given by the probability density function 
Aa " * ” 


W p(W ) = 0.441 v 2 e ^ V ^ . (4.6) 

A a A a 


A plot of Equation 4.6 is given in Figure 4-10. 

Values of k^ , k^, and k^ are given in Figures 4-2, 4-9, and 4-11, 

respectively. The expression on the abscissa of Figure 4-11 is a dimensionless 

quantity. The value of a is 4. 5 if a = x and 7. 5 if a = z or y; A a is the 

a 

distance over which the gust is 50 -percent correlated, and 17 is given in 
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the insert. To complete the definition, a set of discrete gusts of period r is 

selected according to the following criteria: Let r (W. ) be the most probable 

o Act 

period of gust with amplitude W . In Figures 4-5 through 4-8 relationships 
are proposed which relate gust period and gust amplitude. Then 



where n and n . are obtained from known or assumed response character- 
max mm 

istics of the specific wind turbine. The effects of these discrete gusts on a given 
wind turbine are then predicted by means of a deterministic aero-structural 
dynamic model. In this way a statistical description of the wind turbine loads is 
obtained from the statistical description of the wind. 

4. 2. 2. 3 Number of Cycles 

The number of cycles in wind speed which exceed the gust amplitude, 
w a«< a a , h) , in a given year’s exposure to winds having a Rayleigh distribution 

A 

is determined from Figures 4-12 and 4-13 by (the annual mean wind speed at 

height h) and Act (the longitudinal, lateral, or vertical spatial extent of the gust, 

i. e. , a gust of size large enough that it is expected to engulf a distance Aa ) . 

Values of k„ . are given in Figure 4-14. 
lOO! 
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Figure 4-1. Dimensionless turbulence kinetic energy spectra. 
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Figure 4-2. Turbulence intensity <r /W, „ „ (dimensionless), 

w h=10 m ' ' 
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Figure 4-3. Period of gust 50-percent coherent over the horizontal 

distance, A a? = Ax. 
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PERIOD OF GUST, t (s) 



SPATIAL DISTANCE, Aa (m) 


Figure 4-4. Period of gust 50-percent coherent over the lateral 
or vertical distance, Aa = Ay = Az. 
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(DIMENSIONLESS) 
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■I 



HEIGHT, h (m) 


Figure 4-5. Factor for computing the discrete gust amplitude, 

W. = kkW , where k is given in Figure 4-9 

1 5 6 h=10 m 6 & 

and W, „ _ = 5 m s -1 . 

h=10 m 
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(DIMENSIONLESS) 


in 



300 


Figure 4-6. Factor for computing the discrete gust amplitude, 


W. = k k W 


5 6 h=10 m 


, where k is given in Figure 4-9 


and W, = 10 m s 1 . 

h=10 m 
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(DIMENSIONLESS) 



5 10 BO 100 300 


HEIGHT, h (m) 

Figure 4-7. Factor for computing the discrete gust amplitude, 

W. = k_k„W, , , where k„ is given in Figure 4-9 

l 5 6 h=10 m 6 & 

and W, = 15 m s -1 . 

h=10 m 
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(DIMENSIONLESS) 


in 



Figure 4-8. Factor for computing the discrete gust amplitude, 

W. = k k W , where k„ is given in Figure 4-9 

1 5 6 h=lO m 6° & 

andW, > 25 ms" 1 . 

h=10 m 
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Figure 4-9. Factor for adjusting W, „ to W n . 

J b h=10 m h 
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Figure 4-10. Probability density function of the discrete gust amplitude, W 
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Figure 4-11. Factor for estimating effective rms value, for cyclic gusts. 



(DIMENSIONLESS) 



Figure 4-12. Factor for calculating the number of times per year 

the wind speed exceeds the value W . 

A a 
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(DIMENSIONLESS) 


a a Aa(7 W /h 


o> 



a a Aa(7 W /h 

F igure 4-13. Zero crossing factor. 
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SURFACE ROUGHNESS LENGTH, z Q (m) 


Figure 4-14. Scaling factor. 


Chapter 5.0 Wind Direction 


5. 1 Wind Direction Probability 

The cumulative probability of the horizontal wind vector lying within an 
angle between 0 and 0 is given in Figures 5-1 through 5-5 as a function of height, 
h, and surface roughness length, z q (z q values are given in Table 3-1). The 

angle, 0=0, is the direction in which the mean horizontal wind speed is blowing. 
The cumulative probability curves are symmetric such that the probability of 0 
lying within the arc 20 is twice the probability of the angle lying between 0 and 0. 

5. 2 Wind Direction Fluctuations 

The number of times per unit time the wind direction fluctuations exceed 

the angle, 0, is given in Figures 5-6 through 5-9 as a function of height, h, 

surface roughness length, z , and the mean wind speed at a height of 10 m, 

W. in . That is, ° 

h=10 m 


N(0) = k n 


< W h=10 J 


0. 17 


where W, is in meters per second and N(0) is in seconds. Doubling the 

h=10 m v ' 

value of N(0) provides the number of times per unit time that the horizontal wind 

vector fluctuates outside a given arc of 20 (i. e. , the number of times per unit 

time the angle exceeds ±0). These results are relative to a 1-h averaging period. 
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9 (deg) 







Figure 5-3. Cumulative probability distribution of angular 
displacement of mean wind, h = 50 m. 
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Figure 5-4. Cumulative probability distribution of angular 
displacement of mean wind, h = 100 m. 
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Figure 5-5. Cumulative probability distribution of angular 
displacement of mean wind, h = 150 m. 
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(DIMENSIONLESS) 



Figure 5-6. Factor for computing the number of times per unit time the wind 
vector fluctuation exceeds the angle 0 measured from the direction 
of the mean wind vector, z Q = 0.001 m. 
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(DIMENSIONLESS) 



9 (deg) 


Figure 5-7. Factor for computing the number of times per unit time the wind 
vector fluctuation exceeds the angle 0 measured from the direction 
of the mean wind vector, z Q = 0.01 m. 
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(DIMENSIONLESS) 



Figure 5-8. Factor for computing the number of times per unit time the wind 
vector fluctuation exceeds the angle 9 measured from the direction 
of the mean wind vector, z Q = 0.1 m. 
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(DIMENSIONLESS) 



0 (deg) 


Figure 5-9. Factor for computing the number of times per unit time the wind 
vector fluctuation exceeds the angle e measured from the direction 
of the mean wind vector, z Q = 1.0 m. 
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